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Abstract 

In this paper, we establish lower and upper Gaussian bounds for the probability density 
of the mild solution to the stochastic heat equation with multiplicative noise and in any 
space dimension. The driving perturbation is a Gaussian noise which is white in time with 
some spatially homogeneous covariance. These estimates are obtained using tools of the 
Malliavin calculus. The most challenging part is the lower bound, which is obtained by 
adapting a general method developed by Kohatsu-Higa to the underlying spatially homo- 
geneous Gaussian setting. Both lower and upper estimates have the same form: a Gaussian 
density with a variance which is equal to that of the mild solution of the corresponding 
linear equation with additive noise. 
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1 Introduction and main result 



In this paper, we aim to establish Gaussian lower and upper estimates for the probability 
density of the solution to the following stochastic heat equation in R*^: 

^{t,x) - Au{t,x) = b{u{t,x)) +(7{u{t,x))W{t,x), (t, x) G [0, T] X R'^, (1.1) 

with initial condition m(0, x) = uo{x), x € R*^. Here, T > stands for a fixed time horizon, the 
coefficients a, 6 : R — t- R are smooth functions and uq : R'' i— >■ R is assumed to be measurable 
and bounded. As far as the driving perturbation is concerned, we will assume that W{t, x) is 
a Gaussian noise which is white in time and has a spatially homogeneous covariance. This can 
be formally written as: 

E [W{t, x)W{s, y)] = 6{t - s)A{x - y), (1.2) 

where S denotes the Dirac delta function at zero and A is some tempered distribution on R*^ 
which is the Fourier transform of a non-negative tempered measure fi on R'^ (the rigorous 
definition of this Gaussian noise will be given in Section 12. ip . The measure /i is usually called 
the spectral measure of the noise W. 

The solution to equation (jl.ip will be understood in the mild sense, as follows. Let {^t)t>o 
denote the filtration generated by the spatially homogeneous noise W (see again Section 12.11 
for its precise definition). We say that an ^t-adapted process {u{t,x), {t,x) e [0,T] X R"^} 
solves (|l.ip if it satisfies: 



u{t,x) = {T{t)*uo){x)+ / r{t-s,x-y)a{u{s,y))W{ds,dy) 

Jo^Jm^ (1.3) 

+ / / ^{t - s,x - y)b{u{s,y))dyds, 
Jo Jr'' 

where * is the standard convolution product in R*^, and F denotes the fundamental solution 
associated to the heat equation on R*^, that is, the Gaussian kernel of variance 2t: T{t,x) = 

(47ri)~^ exp ■'^^^ , for (t, x) G R+ x R'^. Note that the stochastic integral on the right-hand 

side of ()1.3p can be understood either in the sense of Walsh [30], or using the further extension 
of Dalang [4j (see also [SUE] for another equivalent approach). Indeed, 01 Theorem 13] and 
[71 Theorem 4.3] imply that equation ()1.3p has a unique solution which is L^-continuous and 
satisfies, for all p > 1: 

sup E [\u{t,x)\P] < +0O. 

(t,a')G[0,T]xR'* 

Let us point out that, in the above-mentioned results, the fundamental solution F and the 
noise W are related as follows: 

$(r) := r [ \^T{tm\^fi{dOdt < +00. (1.4) 
Jo Jr'' 

This quantity measures the variance of the stochastic integral in (|1.3p (indeed, it is the variance 
itself when £7 = 1), therefore it is natural that $(t) will play an important role in the Gaussian 
lower and upper bounds for the density of the random variable u{t,x). Moreover, it has been 
proved in ^ Example 2] that condition (|1.4p is satisfied if and only if: 

^--pi^ /i(de) < +00. (1.5) 
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We also remark that the stochastic heat equation (jl.ip has also been studied in the more 
abstract framework of Da Prato and Zabczyk [B] and, in this sense, we refer the reader to 
|24j and references therein. Nevertheless, in the case of our spatially homogeneous noise, the 
solution in that more abstract setting could be obtained from the solution to equation (|1.3p 
(see O Sec. 4.5]). 

The techniques of the Malliavin calculus have been applied to equation (jl.Sp in the papers 
[IHET]. Precisely, [211 Theorem 6.2] states that, if the coefficients b and a are ^°°-functions 
with bounded derivatives of order greater than or equal to one, the diffusion coefficient is 
non-degenerate (i.e. |o"(z)| > c > for all z G M), and ()1.5p is satisfied, then for each {t,x) € 
(0,r] X W^, the random variable u{t,x) has a density pt^x (see also Theorem 13.21 below). 
Moreover, in the recent paper [20], the strict positivity of this density has been established 
under a "^^-condition on the density and the additional condition of a being bounded. 

Our aim in this paper is to go a step further and prove the following theorem: 

Theorem 1.1. Assume that condition (jl.Sp is satisfied and a,b ^fj°°(M) bounded and 

bounded derivatives). Moreover, suppose that \o'{z)\ > c > 0, for all z € M. Then, for every 
{t, x) € (0, T] X W^, the law of the random variable u{t, x) has a density pt^x satisfying, for 
all y G M; 

where Fq = {T(t) *uo){x) and ci, C2, C3, Ci, C2 are positive constants that only depend on T, a 
and b. 

One of the interests of these type of bounds is to understand the behavior of the density 
when y is large and t is small. In both cases, one obtains the same upper and lower behavior for 
the density, that is, a Gaussian density with a variance which is equal to that of the stochastic 
integral term in the mild form of the linear equation. We observe that this variance does not 
depend on x due to the spatially homogeneous structure of the noise. 

In order to prove our main result, we will apply the techniques of the Malliavin calculus, 
for which we refer the reader to [18] and [28]. Obtaining lower and upper Gaussian bounds 
for solutions to non-linear stochastic equations using the Malliavin calculus has been a current 
subject of research in the last twenty year. Precisely, the expression for the density arising from 
the integration-by-parts formula of the Malliavin calculus provides a direct way for obtaining 
an upper Gaussian-type bound for the density. Indeed, ones applies Holder's inequality, and 
then combines the exponential martingale inequality together with estimates for the Malliavin 
norms of the derivative and the Malliavin matrix. This is a well-known method that has been 
applied in many situations (see for instance [Hi |5] ) . We will also apply this technique to the 
density of our stochastic heat equation in order to show the upper bound in Theorem 11.11 (see 
Section 5). 

On the other hand, to obtain Gaussian lower bounds for some classes of Wiener functionals 
turns out to be a more difficult and challenging issue. In this sense, the pioneering work is 
the article by Kusuoka and Stroock [T2|, where the techniques of the Malliavin calculus have 
been applied to obtain a Gaussian lower estimate for the density of a uniformly hypoelliptic 
diffusion whose drift is a smooth combination of its diffusion coefficient. Later on, in |lUj . 
Kohatsu-Higa took some of Kusuoka and Stroock's ideas and constructed a general method 
to prove that the density of a multidimensional functional of the Wiener sheet in [0, T] x 
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admits a Gaussian-type lower bound. Then, still in [TO], the author applies his method to a 
one-dimensional stochastic heat equation in [0, 1] driven by the space-time white noise, and 
obtains a lower estimate for the density of the form: 

1 ( \y-Fo\' \ 
Cit 4 exp — , 

V C2t-2 J 

where Fq denotes the contribution of the initial condition. This is the bound we would get if in 
equation (jl.ip we let d = 1 and W be the space-time white noise. Indeed, this case corresponds 
to take A = (5 in ()1.2p . therefore the spectral measure is the Lebesgue measure on and 
^{t) = C t2 . As we will explain below, in the present paper we will use Kohatsu-Higa's method 
adapted to our spatially homogeneous Gaussian setting. The same author applied his method 
in [11] to obtain Gaussian lower bounds for the density of uniformly elliptic non-homogeneous 
diffusions. Another important case to which the method of [TO] has been applied corresponds 
to a two-dimensional diffusion, which is equivalent to deal with a reduced stochastic wave 
equation in spatial dimension one, a problem which has been tackled in |6j. Moreover, let 
us also mention that the ideas of [TO] have been further developed by Bally in [1] in order 
to deal with more general diffusion processes, namely locally elliptic Ito processes, and this 
has been applied for instance in [9j. Eventually, in [H] Bally and Kohatsu-Higa have recently 
combined their ideas in order to obtain lower bounds for the density a class of hypoelliptic 
two-dimensional diffusions, with some applications to mathematical finance. 

The increasing interest in finding Gaussian lower estimates for Wiener functionals has 
produced three very recent new approaches, all based again on Malliavin calculus techniques. 
First, in [T7j the authors provide sufficient conditions on a random variable in the Wiener 
space such that its density exists an admits an explicit formula, from which one can study 
possible Gaussian lower and upper bounds. This result has been applied in |22l [23] to our 
stochastic heat equation (jl.ip in the case where a = 1. Precisely, Theorem 1 and Example 
8] imply that, if b is of class with bounded derivative and condition (jl.Sp is fulfilled, then, 
for sufficiently small t, x) has a density pt^x satisfying, for almost all z E M: 

where Mt^x = E(u{t,x)) (see [221 Theorem 4.4] for a similar result which is valid for all t 
but is not optimal). Compared to Theorem 11.11 one the one hand, we point out that our 
result is valid for a general a, arbitrary time T > and our estimates look somehow more 
Gaussian. On the other hand, the general method that we present in Section 12.21 requires the 
underlying random variable to be smooth in the Malliavin sense, and this forces to consider 
a smooth coefficient b. We also remark that, even though the results of [2.3] are also valid for 
a more general class of SPDEs with additive noise (such as the stochastic wave equation in 
space dimension d G {1,2,3}), Nourdin and Viens' method does not seem to be suitable for 
multiplicative noise settings. 

A second recent method for deriving Gaussian-type lower estimates for multidimensional 
Wiener functionals has been obtained by Malliavin and E.Nualart in [13] (see [TO] for the one- 
dimensional counterpart). This technique is based on an exponential moment condition on the 
divergence of a covering vector field associated to the underlying Wiener functional, and has 
been applied in ^19j to a one-dimensional diffusion. 
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Last, but not least, in the recent paper [2], Bally and Caramellino develop another method 
to obtain lower bounds for multidimensional Wiener functionals based on the Riesz transform. 

As we have already mentioned before, in the present paper we will apply the methodology 
of Kohatsu-Higa [TU]. For this, first we will need to extend the general result [lOj Theorem 
5] on Gaussian lower bounds for uniformly elliptic random vectors from the space-time white 
noise framework to the case of functionals of our Gaussian spatially homogeneous noise (see 
Theorem \2.3\i . This will be done in Section [21 after having precisely described the Gaussian 
setting which we will work in. The extension to such a general case turns out to be quite 
straightforward, since it essentially requires to replace norms in L^([0,T] x A), with A C R'^, 
by those in L'^{[0,T]; J^), where is the Hilbert space that can be naturally associated to 
the spatial structure of the noise (see Section 12.11 below) . 

In Section [3l we will recall the main results on differentiability in the Malliavin sense 
and existence and smoothness of the density applied to our stochastic heat equation (|1.3p . 
Moreover, we will prove a technical and useful result which provides a uniform estimate for the 
conditional norm of the iterated Malliavin derivative of the solution on a small time interval. 

Section m is devoted to apply the general result Theorem 12. 3fo f Section [212] to the stochastic 
heat equation ()1.3p . to end up with the lower bound in Theorem 11.11 That is, one needs to 
show that the solution u{t, x) defines a uniformly elliptic random variable in the sense of 
Definition 12.21 Although the proof's structure is similar to that of |lQi Theorem 10], the 
analysis in our setting becomes much more involved because of the spatial covariance structure 
of the underlying Wiener noise. As mentioned before, the upper bound in Theorem 1 1 . 1 1 will be 
proved in Section [5l 

Eventually, we have also added an appendix where, first, we recall some facts concerning 
Hilbert-space- valued stochastic and pathwise integrals and their conditional moment estimates, 
and, secondly, we state and prove a technical result which has been applied in Section [H 

As usual, we shall denote by c, C any positive constants whose dependence will be clear 
from the context and their values may change from one line to another. 

2 General theory on lower bounds for densities 

This section is devoted to extend Kohatsu-Higa's result Theorem 5] on lower bounds for 
the density of a uniformly elliptic random vector to a more general Gaussian space, namely 
the one determined by a Gaussian random noise on [0, T] x M'^ which is white in time and has 
a non-trivial homogeneous structure in space. For this, first we will rigorously introduce the 
Gaussian noise and the Malliavin calculus framework associated to it and needed in the sequel. 

2.1 Gaussian context and Malliavin calculus 

Our spatially homogeneous Gaussian noise is described as follows. On a complete probability 
space (r2,^,P), let W = {W{ip),ip G '^q°°(M+ x R^)} be a zero mean Gaussian family of 
random variables indexed by functions with compact support with covariance functional 
given by 

E[W(ip)Wii;)]= [ dt [ A{dx)(ip{t,i^)*i;{t,i.)){x), ip,'tp e '^^{R+ x R"^). (2.1) 

Jo JRd V / 

Here, A denotes a non-negative and non-negative definite tempered measure on R'^, * stands for 
the convolution product, the symbol ★ denotes the spatial variable and i/j(t, x) := ip(t, —x). For 
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such a Gaussian process to exist, it is necessary and sufficient that the covariance functional 
is non-negative definite and this is equivalent to the fact that A is the Fourier transform of a 
non- negative tempered measure /i on M.'^ (see |29^ Chap. VII, Theoreme XVII]). The measure 
/i is usually called the spectral measure of the noise W. By definition of the Fourier transform 
of tempered distributions, A = JF^ means that, for all <p belonging to the space ^(K'^) of 
rapidly decreasing functions, 

cp{x)A{dx) = [ ,^mm)- 



< +00. 



Moreover, for some integer m > 1 it holds that 



I 



(1 + 



\2\m 



Elementary properties of the convolution and Fourier transform show that covariance ()2.ip can 
be written in terms of the measure fi, as follows: 



In particular, we obtain that 



/'OO f' 

Jo JR<< 



/o Jr^ 

Example 2.1. Assume that the measure A is absolutely continuous with respect to Lebesgue 
measure on M*^ with density f. Then, the covariance functional (|2.ip reads 



00 



dt dx I dyip{t,x)f{x-y)tp{t,y), 





which clearly exhibits the spatially homogeneous nature of the noise. The space-time white 
noise would correspond to the case where f is the Dirac delta at the origin. 

We note that the above-defined kind of noise has been widely used as a random perturbation 
for several classes of SPDEs (see for instance [15l HI [Ml l2H] ) • 

At this point, we can describe the Gaussian framework which is naturally associated to our 
noise W . Precisely, let be the completion of the Schwartz space 5^{W^) endowed with the 
semi-inner product 



{(^iA2).^= {(t>i * (t>2){x) A{dx) = ^MO^MOKdO, 4>iA2ey{R''). 

Notice that the Hilbert space may contain distributions (see [H Example 6]). Fix T > and 
define J^t = -^^([0,7"];^). Then, the family W can be extended to J^t, so that we end up 
with a family of centered Gaussian random variables, still denoted by = {W{g),g € J^}, 
satisfying that E [1^(511) 1^(5(2)] = {gi,g2),^T^ 91-, 92 £ -J^T (see for instance [3 Lemma 

2.4] and the explanation thereafter). 

The family W defines as isonormal Gaussian process on the Hilbert space and we shall 
use the differential Malliavin calculus based on it (see, for instance, [ISIEB])- We denote the 
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Malliavin derivative by D, which is a closed and unbounded operator defined in L^(r2) and 
taking values in L'^(Q]J^), whose domain is denoted by B-*^'^. More general, for any m > 1, 
the domain of the iterated derivative D"^ in LP(0) is denoted by D^'^, for any p > 2, and we 
recah that D"^ takes values in LP(0;^®™). As usual, we set = np>i nm>i B"^'P. The 
space ~™ " 



can also be seen as the completion of the set of smooth functionals with respect to 



the semi-norm 



m 

\F\\m,p = {E[\F\P] + J2^ 



ID^FI 



®3 



D"^F{r) is an 



For any differentiable random variable F and any r 
element of Jf®™ which will be denoted by L>™F. 

We define the Malliavin matrix of a fc-dimensional random vector F € { 
{{DFi, DFj) j^^)Kij<k- We will say that a A:-dimensional random vector F is smooth if each 
of its components belongs to B°°, and we will say that a smooth random vector F is non- 
degenerate if (det 7f)~^ € r\p>iL^{Q). It is well-known that a non-degenerate random vector 
has a density (cf. [ig Theorem 2.1.4]). 

Let {'-!^t)t>o denote the cj-field generated by the random variables {Ws{h), h € J^, < s < 
t} and the P-null sets, where Wt{h) := W{l[Qt]h), for any t > 0, h G J^. Notice that this 
family defines a standard cylindrical Wiener process on the Hilbert space J^. We define the 
predictable u-field as the cr-field in $7 x [0,T] generated by the sets {{s,t] x A,0 < s < t < 

T,Ae^s}. 

As in |10[ Section 2], one can define the conditional versions of the above Malliavin norms 
and spaces (see also pTOl For all < s < t < T, we set Jifs,t = -^^([s, t]; For any 
integer m > and p > 1, we define the seminorm: 



I llm,p 



m 

{E,[|Fn + ^E, 

.7 = 1 



ID^FI 



where Es[-] = E[-|^5]. We will also write Ps{-} = P{-|=^s}. Completing the space of smooth 
functionals with respect to this seminorm we obtain the space B^^^. We write L^j(Q) for 
B°f . We say that F G W^f if F E B™f and ||F|[^*p G ng>iL9(0). If we let now F be a 
fc-dimensional random vector, we define its associated conditional Malliavin matrix as 7^* = 
{{DFi, DFj) t)i<i,j<k- We say that F = {Fi, . . . , F^) is conditionally non-degenerate in [s, t] 
if Fi e 



r\ni>o,p>i'^s,t 1 for alH = 1, . . . , /c, and (det 7^ ) ^ G np>iL^g^^{Q). 



2.2 The general result 

In order to state the main result of this section, we need to define what we understand by a 
uniformly elliptic random vector in our context. 

Definition 2.2. Let F be a smooth ^t-in^osurahle k- dimensional random vector. We say that 
F is uniformly elliptic if there exists an element g G such that ||5'(s)||jf' > for almost all 
s, and an € > such that, for any sequence of partitions tin = {Q = t^ < ti < ■ ■ ■ < t^ = T} 
whose norm is smaller than e and WttnW = sup{i„ — n = 1, A^} converges to zero as 

N —7- 00, the following conditions are satisfied: 
Define 

0<A„_i(5):= / ||c/(s)|[^ds < 00, n = l,...,A^. 

''tn-l 
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(HI) (Approximation property) There exists a sequence of smooth random vectors Fq, Fi,...,Fn 

such that F/v = F, Fn is ^t„-'nT'Casurable and, for any n = 1, . . . , N , Fn is conditionally 

non- degenerate in [tn-ij^n]- Moreover, for some 7 > and for each F^ and K ^ 'H, 

— — K 

there exists a random, vector Fn = F^ of the form 

Fn = A„_i(5)(^+i)^Z„ + Fn-l + I^ih) + G^, 

where the various terms in this formula and the relationship with F^ are as follows. 

{Zn, n = 1, . . . , N) is an i.i.d. sequence of k-dimensional N{0, Id) random vectors which 
are independent of the Gaussian family W , where Id denotes the k-dimensional identity 
matrix. The random vectors are J^t^ -measurable and belong to ^t„-i,tn- Eventually, 
I'^(h) denotes a random vector whose components are of the form 

{r{h)Y=r f h'{s,y)W{ds,dy), e = l,...,k, 

where, for each £, h^ is a smooth ^^^ -^ -pre(iicta6Ze processes with values in ^tn-i,tn o-nd, 
for any m G N and p > 1, there exists a constant G such that: 

||-Fn||m,p + sup \\h^\\Ml^_^,tM ^ C"' 

for any I = 1,. . . ,k. 

Furthermore, the following four conditions are satisfied for the approximating sequence 
Fn, for all n = 1, . . . , N: 

(H2a) For any m G N and p > I: 

< CA„_i(5)(^+i)^ a.s. 

(H2b) For any p > 1: 

{E*„_Jdet(7l.:-^'*'^)-^]}^/^ ^ C^n-i{g)-' a.s. 
(H2c) Let A = (a^^q) be the k x k matrix defined by 

a,,, = A„_i(<7)-i / {h'{s), hiis))^ds. 

There exist strictly positive constants G\ and C2 such that, for all ( eM/^, 

CiCC, > C^^C > CsC^C, a.s. 
(H2d) There exists 5 > Q such that, for all m gN and p > 1: 

\\G^\t:p'''''<CAn-i{g)'/'+' a.s. 
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In hypothesis (HI) above, the Gaussian random variable is indeed defined as follows. 
One enlarges the underlying probability space in order to include the increments of another 
cylindrical Wiener process \W t(h\ t >0,h £ J^} independent of W. Then, = Wn+i{h) — 
Wn{h), where h is any element of a complete orthonormal system on Jif. In this case, the 
expanded filtration is given by = V a{Ws{h), s < n + l,tn < t,h £ Jif}. All the norms 
and expectations considered so far have to be understood in the enlarged probability space. 
We remark that, as it has been explained in [10^ p. 436] for the space-time white noise case, 
the introduction of the term An-i{g)^^^^^"' Zn is necessary in order to ensure that F„ has a 
density. 

On the other hand, note that condition (H2c) is the ingredient that most directly reflects 
the uniformly elliptic condition for a random vector on the Wiener space. 

The next theorem establishes a Gaussian lower bound for the probability density of a 
uniformly elliptic random vector. Its proof follows exactly the same steps as that of Theorem 5 
in [TO], where the same result has been proved in a Gaussian setting associated to the Hilbert 
space L2([0,r] X A), where A C M'^ (that is, the space-time white noise). In our case, one 
essentially needs to replace the norms in the latter space by those in Mr = L'^{[0,T]; J^), 
in such a way that the proof follows with minimal changes. Nevertheless, for the sake of 
completeness, we believe that it is worth reminding the main steps of the proof of this theorem, 
for the reader will have a much clearer understanding of the long list of assumptions stated 
above. 

Theorem 2.3. Let F be a k- dimensional uniformly elliptic random vector and denote by pp 
its probability density. Then, there exists a constant M > that depends on all the constants 
of Definition \2.S\ such that: 

PF{y) > M\\g\\J/\^pU^^^-^\ for ally GR", 

\ \\9\\ / 

where Fq is the first element in the sequence of hypothesis (HI). 
Sketch of the proof. The proof is divided in four steps. 

Step 1. First, as we have already mentioned, the fact that in the definition of F^ we include 
the term involving the random variable Z„ let us prove that Fn is conditionally non-degenerate 
in Therefore, it has a smooth conditional density with respect to ^tn-i- 

Next, one proves that assumptions (H2a)-(H2d) imply the following local lower bound 
for this conditional density: there exist constants M, c and rj such that, if lS.n-i{g) < r] and 
yeAn = {yeR'': \\y - F„_i|| < cA^-iig)^/^}, then 

E.„-.[^.(F.)]> ^^^J^^^^,,, a.s., (2.2) 

where 6y denotes the Dirac delta function. This is proved using the expression for the density 
that follows from the integration-by-parts formula of the Malliavin calculus. After normalizing 
this density (i.e. dividing by A„_i 

(5)^/2)^ the Taylor expansion of the delta function around 
the non-degenerate random vector An-i{g)~^^'^ I"' (h) is considered. Condition (H2d) yields 
that the terms of higher order, that is, those that are concentrated in the smooth random 
vector , are of order An-i{g)^ , so that they are negligible with respect to the first term 
I^{h). Using the fact that (H2c) is equivalent to saying that the density of the term I"'{h) 
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corresponds to a non-degenerate Gaussian random variable, the lower bound (j2.2p for the 
conditional density of Fn is obtained. 

Step 2. The second step consists in proving that, if Fn{p) ■= pFn + (1 — p £ [0, 1], then 

for any p > 1 there exists a constant C such that: 

sup {E,„_,[det(7i"Y;;")-^]}^/^ < CA^.,{g)-' a.s. (2.3) 

The proof of this inequality is very technical (see [10, Proposition 12]), but (|2.3p is an important 
ingredient, together with ()2.2p . to end up with a lower bound for the conditional density of Fn 
with respect to ^t„-i- This is explained in the next step. 

Step 3. Next, one shows that (H2a), (j2.2p and (|2.3p imply that there exist positive con- 
stants c, M, a > l,r/, and random variables C„ G ''^ = 0, ...,A^ — 1, satisfying that 
sup„=o Af I^(|C'n|) < Af, such that if /S.n-i{g) < rj and y G ^n, then 

Et„.,[<5,(F„)] > - C„-i(a;)A„_i(5)", (2.4) 

for almost all w G fi. This lower bound is proved by writing 

Ei„_J<5,(F„)] > _^-l^ + Ei„_J<5,(F„)-5,(F„)], 

and then finding an upper bound of the second term on the right-hand side using (H2a) and 

(123]). 



Step 4- Finally, one concludes the desired lower bound of Theorem 12.31 using ()2.4p as in [10 
Theorem 2]. The main idea is to use Fubini's theorem to write 



nMF)] = ■■■ nSy{F)5y^_,iFN-i) ■ ■ ■ 5y,iFi)]dyi ■ ■ ■ dvN- 



1- 



Then one iteratively applies the lower bound ()2.4p using conditional densities, and finally use 
a localization procedure to conclude the proof. □ 

3 The stochastic heat equation 

In this section, we will recall some known facts about the stochastic heat equation on which 
will be needed in the sequel. We will also prove an estimate involving the iterated Malliavin 
derivative of the solution which, as far as we know, does not seem to exist in the literature (see 
Lemma 13.31 below) . 

We remind that the mild solution to the stochastic heat equation p.ip is given by the 
=^t-adapted process {u{t,x), {t,x) e [0,T] X M'^} that satisfies: 

u(t,x) = (r(t)*uo)(x)+ / / T{t-s,x-y)a{u{s,y))W{ds,dy) 

JoU^ (3.1) 

+ 11 ^{t - s,x - y)b{u{s,y))dyds, 
Jo JR<* 
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where T{t) is the Gaussian kernel with variance 2t and the following condition is fulfilled: 

$(r) = r [ \^rm)\^ fiid^dt < oc. 

Jo Jw^ 

As mentioned in the Introduction, this is equivalent to say that /jgd(l + < co. 

Moreover, similar to [14^ Lemma 3.1], one easily proves that, for all < ri < r2 < T: 

C{T2-n)<r[ \^Tm)\^KdOdt, (3.2) 

for some positive constant C depending on T. 

In order to apply the techniques of the Malliavin calculus to the solution of p.3|) . let us 
consider the Gaussian context described in Section [2m That is, let {W{h), h € Mr} be the 
isonormal Gaussian process on the Hilbert space Mr = L'^{[0,T];J^) defined therein. Then, 
the following result is a direct consequence of |14^ Proposition 2.4], [27, Theorem 1] and |2H 
Proposition 6.1]. For the statement, we will use the following notation: for any m G N, set 
s := {si,...,Sm) G [0,r]"", z := G {R'^)"^, s{j) := {si, . . . , sj-i, sj+i, . . . , Sm) 

(resp. z{j)), and, for any function / and variable X for which it makes sense, set 

A™(/,X) := D"'f{X) - /'(X)D™X 

Note that A"^{f, X) = for m = 1 and, if m > 1, it only involves iterated Malliavin derivatives 
up to order m — 1. 

Proposition 3.1. Assume that (|1.5|) is satisfied and a,b & 'i^'°°(M) and their derivatives of 
order greater than or equal to one are bounded. Then, for all (i, x) G [0, T] x R'^, the random 
variable u{t,x) belongs to 0°° . Furthermore, for any m G N and p> 1, the iterated Malliavin 
derivative D"^u{t,x) satisfies the following equation in {Q.; .^^^) : 

D"'u{t,x) = Z"^(t,x) 

+ 1^ r{t -s,x- y)[A"^(a, u{s, y)) + Z)™n(s, y)a'{u{s, y))]W{ds, dy) ^3 3^ 

+ /* / T{t-s,x-y)[JV^{b,u{s,y)) + D'^u{s,y)b'{u{s,y))\dyds, 
Jo JR'* 

where Z"^(t,x) is the element of {Q; J^r^"^) given by 

m 

Z^{t,x)s,2 = J2^{t- Sj,x- zj)D^^~].^^^a{u{sj,z,)). 
i=i 

We remark that the Hilbert-space- valued stochastic and pathwise integrals in equation ()3.3p 
are understood as it has been described in Section lA.ll 

As far as the existence of a smooth density is concerned, we have the following result (see 
[211 Theorem 6.2]): 

Theorem 3.2. Assume that (|1.5|) is satisfied and o", 6 G 'rf°°{M) and their derivatives of order 
greater than or equal to one are bounded. Moreover, suppose that \o'{z)\ > c > 0, for all ^ G M. 
Then, for every {t, x) G (0, T] x M'^, the law of the random variable u{t, x) has a density. 
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The following technical result, which will be used in the proof of Theorem II. 1^ exhibits 
an almost sure estimate for the conditional moment of the iterated Malliavin derivative of u 
in a small time interval. As will be explained in Remark 13.41 this result is still valid for a 
slightly more general class of SPDEs, such as the stochastic wave equation in space dimension 
d € {1, 2, 3}. Nevertheless, for the sake of simplicity, we will focus either the statement and its 
proof on our stochastic heat equation (jl.3p . 



Lemma 3.3. Let 0<a<e<T, m€N and p > 1. Assume that the coefficients 6, cr : M — )• M 
belong to ^°°(M) and all their derivatives of order greater than or equal to one are bounded. 
Then, there exists a positive constant C, which is independent of a and e, such that, for all 
6 e (0, e - a] ; 



sup Ea ( ||L»'"u(r,y)ff,55„ 

'e — S,e 



mp 



a.s.. 



(T,y)e[e—5,e]xl 

where we remind that J^^^^ denotes the Hilhert space LF'{[e — S, e]; J^^"^) and, for all t>0, 

ft 



JM.'' 



Proof. We will proceed by induction with respect to m G N. First, let us observe that the case 
m = 1 has been proved in [U Lemma 2.5] (see also [271 Lemma 5]). Suppose now that the 
statement holds for any j = 1, . . . ,m — 1, and let us check its veracity for j = m. 

Let e — 6 < t < e and x G M'^. Then, the conditioned norm of the Malliavin derivative 
W^ui^t, x) can be decomposed as follows: 



Ea \\D"'uit,x) 



|2p 



<CiBi + B2 + B3), a.s. 



with 



Y^r{t-sj,x-*)Df-^a{uis„*)) 



ds 



(here, if we formally denote by (zi, . . . , Zm) the variables of J^f®'^, the symbol ★ corresponds 

to Zj), 



B2 = Ea 



B:> = Ea 



t 

Jm.'^ 
t 



T{t -s,x- y)[A™(a, n(s, y)) - D^u{s, y)a'{u{s, y))]W{ds, dy) 

2p 



2p 



e — (5,6/ 



JK.'^ 



T{t -s,x- y)[A"(6, u{s, y)) - D^u{s, y)b'{u{s, y))] dyds 



e — 5,e/ 



Let us start with the study of the term Bi. First, note that we must have that e — 5 < Sj < t, 
thus 



Bi<CY,Ea / 

i=i v-^^- 



ds{j) T{t - Sj,x - *)L»^^-)V(^/(s,-,*)) 



2 
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At this point, we can proceed as in the proof of [27', Lemma 2] (see p. 173 therein), so that we 
can infer that 



Bi<c(r J{t-r)dry f snvEa{\\D^-^a{u{r,y))f^ 



1) J{t-r)dr, 



where we have used the notation 



J(r)= / |^r(r)(e)|V(de). 

Precisely, we have used the fact that F is a smooth function, and then apphed Holder's and 
Cauchy-Schwarz inequalities. Hence, we have that 



Bi <C{mr sup f ||I)"^-VKr,y))f^^,_„) 



(3.4) 



In order to bound the above supremum, one applies the Leibniz rule for the iterated Malliavin 
derivative, the smoothness assumptions on a, Holder's inequality and the induction hypothesis, 
altogether yielding 



sup Ea (\\D^-'cT{u{r,y))fP < C(cI>(^))('"-l)^ 



a.s. 



Plugging this bound in ()3.4p . we end up with 

Bi < C($(5))™P, a.s. 



(3.5) 



Next, we will deal with the term B2, which will be essentially bounded by means of Lemma 
lA.H as follows: 



B2<C (^J^ J{t - r) dr^ 



p-i 



le-5 



sup Ea ( |IA"^(^7,7x(s,y))||5^„ ) + sup Ea ( \\D^u{s,y)\ 



|j2p 

II jt^®m 



J{t — s)ds. 



Owing again to the the Leibniz rule for the Malliavin derivative and noting that A™ only 
involves Malliavin derivatives up to order m — 1, one makes use of the induction hypothesis to 
infer that 



sup E, (|IA-(a,n(.,y))|| J«„ ) < C(cl>(5))™^ 



a.s. 



Hence, 



B2 < cmryr-' f imr" + sup e^ f p™^.(s,y)f^«„, \ 



< Ci f {mr" + sup Ea (\\D^u{t, 

Je~S [ (T,j/)e[e-5,s]xR<* V 

almost surely, where Ci denotes some positive constant. 



||2p 

e — S,e 



J{t — s)ds 
J{t — s)ds^ 



(3.6) 
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Eventually, using similar arguments, we can show that the term is bounded above by: 



C 



ie~5 



T{t - s,x-y) 



En 



|A-(6,n(.,y))|| 



2p 



+ Ea \\D^u{s, 



|2p 

1^5 



dyds 



< C 



imr" + sup Ea (\\D^u{T,y)\\^^,^,„ ) 

{T,y)£[e-S,s]xR'i \ '^e-S^e J 



ds, a.s. 



(3.7) 



Here, we have also used that r(s, y) dy is uniformly bounded with respect to s. 
Set 

F{t):= sup Ea(\\D^u{s,y)\\^p), te[e-6,e]. 

{s,y)e[e~S,t]xR'i \ -^e-S^e J 

Then, ([33])- ([321) imply that 

F{t) < C2md)y''f + Ci f [{^6)rP + F{s)]{J{t -s) + l)ds 

Je-5 



a.s., 



where Ci and C2 are some positive constants. We conclude the proof by applying Gronwall's 
lemma ^ Lemma 15]. □ 

Remark 3.4. Lemma 13.31 still remains valid for a more general class of SPDEs, namely for 
those that have been considered in the paper [4] (see also [21] )• In these references, an SPDE 
driven by a linear second-order partial differential operator has been considered, where one 
assumes that the corresponding fundamental solution F satisfies the following: for all s, T{s) 
is a non-negative measure which defines a distribution with rapid decrease such that condition 
(fO|) is fulfilled and 

sup r(s,M^) < +00. 

0<s<T 

As explained in [H Section 3] , together with the stochastic heat equation, the stochastic wave 
equation in space dimension d € {1,2,3} is another example of such a type of equations. 
Finally, we point out that the proof of Lemma 13.31 in such a general setting would require a 
smoothing procedure of T in terms of an approximation of the identity, which makes the proof 
slightly longer and more technical; this argument has been used for instance in \27\ Lemma 5]. 



4 Proof of the lower bound 

In this section, we prove the lower bound in the statement of Theorem 11.11 For this, we are 
going to show that u{t, x) is a uniformly elliptic random variable in the sense of Definition 12.21 
Then, an application of Theorem 12.31 will give us the desired lower bound. 

To begin with, we fix (t, x) G (0, T] x M.'^, we consider a partition = to < ti < ■ ■ ■ < t]\f = t 
whose norm converges to zero, and define: 



{r{t)*uo){x)+ I I r{t-s,x-y)a{u{s,y))W{ds,dy) 

r(t - s,x - y)h{u{s, y))dsdy. 
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Clearly, F„ is -measurable, for all n = 0, and note that Fq = {T{t)*UQ){x). Moreover, 
Fn belongs to and, for all m G N and p > 1, the norm ||-Fn||m,p can be uniformly bounded 
with respect to {t,x) £ (0,T] x R'^ (see [111 Proposition 2.4], [271 Theorem 1], and also |2ll 
Proposition 6.1]). 

The local variance of the random variable u{t, x) will be measured through the function 
g{s) := r(t — s). Then, observe that we have: 



A„-i(<?) 



tn-l 



and this quantity is clearly positive for all n (see ()1.5p ). 

We next prove that condition (H2b) is satisfied, which in particular implies that Fn is 
conditionally non-degenerate in [t„_i,t„]. Recall that we are assuming that the coefficients b 
and a belong to 'r^°°(M). 

Proposition 4.1. For any p > 0, there exists a constant C such that: 



DrFn\\,%-'dr 



tn-1 



<CA„_i(5)-f a.s. 



(4.1) 



Proof. It follows similarly as the proof of Theorem 6.2 in [21j. More precisely, it suffices to 
show that for any (? > 2, there exists eo = ^o{q) > such that, for all e < eo: 

Pt.-i jA;!i(5) T" \\DrFj%dr < e| < Ce'^ a.s. (4.2) 
Indeed, if we set X = A-^^{g) //"_ WDrFj^.. dr, then we have: 

Choosing q sufficiently large in (j4.2p (namely q > p), we conclude that (|4.ip is fulfilled, and 
hence the statement of (H2b). 

We next prove (|4.2p . For any 6 G (0,t„ — tn-i), we have that 



tn-l 



\DrFj^ dr > -^lo - (2/i + 2/2), 



where 



h 
h 



tn~5 
in 

tn—S 
in 

in—S 



|^r(t-r)(OlV(de)d^, 



r(t — s,x — z)a'{u{s, z))Dru{s, z) W{ds, dz) 



dr, 



T{t — s,x — z)h'{u{s, z))Dru{s, z) dzds 



dr, 



and where we have denoted by ctq the constant such that \(y{v)\ > do > 0, for all u G M. 
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We next bound the pth moments of Ii and I2 for p > 1. Owing to Holder's inequality and 
Lemmas lA.ll and \'6.'6\ we get: 



tn~S 



tn—S 



r(t — s,x — z) 



X a'{u{s, z))Dru{s, z) W{ds, dz) 



tn—5 



2p ■ 

dr 

-if - 

r(t — s,x — z) 



X cr'(n(s, z))Dt„-ru{s, z) W{ds, dz) 



2p 



dr 



< 5^-'ll sup Ei_, 

se[o,5] 



\Dt„-.u{tn - S,z) 



|2p 



Similarly, appealing to the (conditional) Holder's inequality, and Lemmas IA.2I and 13.31 we 
have: 



Et„_J|/2r]<5^-^E,„_, 



tn—5 



T{t — s,x — z)b'{u{s, z))Dt„^ru{s, z) dzds 



2p 



dr 



< 5P-^ 



-1 



tn-S 



\ 2p 



r{t — s, z) dzds j sup ^tn-i 



\Dt„-.u{tn - S,z) 



|2p 



\ 2p 

r{t- s,z)dzds] ^{6)P. 
J 



Observe that, because L is a Gaussian density. 



tn—5 



T{t — s, z) dzds < C 5. 



(4.3) 



Putting together the bounds obtained so far, we have proved that: 

Thus, the (conditional) Chebyschev's inequality yields 

Pt„_jA-\{g) \\DrFn\\%.dr<e 

At this point, we choose 5 = 5{e) such that ^ \(g) ~ Thus, condition (|3.2p implies that 
> ^ '^i(g) ' t^^t is (5 < Ce. Hence, plugging this bound in (|4.3p . we obtain: 

Y>t^K\{g) r \\DrF^t^dr<e]<CI^^eP-^{ll + ll) 

< CeP-\ 
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In order to obtain (j4.2p . it suffices to choose p sufficiently large such that p—1 > q. The proof 
of (|4.ip is now complete. □ 

Remark 4.2. In the above proof, we have used the lower bound ()3.2p . This has prevented 
us from proving our main result Theorem 11.11 for other type of SPDEs, such as the stochastic 
wave equation (see Remark 13. 4p . Indeed, for the latter SPDE, we do not have a kind of time 
homogeneous lower bound of the form (|3.2p . which has been a key point in order to conclude 
the proof of Proposition 14.11 

In order to define the approximation sequence Fn, we proceed similarly as in [10\ p. 442]. 
Precisely, we observe first that 

Fn-Fn-i= r [ r{t-s,x-y)a{u{s,y))W{ds,dy) 

+ / / T{t-s,x-y)b{u{s,y))dsdy. 

Our aim is to find a Taylor-type expansion of the two terms above. This will be done by 
applying the mean value theorem to the functions a and b around the point m„_i(s, y) defined 
by 

n„_.i(s,y)= / T{s,y - z)uo{z)dz + f f T{s - r,y - z)a{u{r, z))W{dr,dz) 

jR'i JO jR'i 

rtn-i r 

+ / T{s — r,y — z)b{u{r, z)) dzdr, 

Jo JR'i 

where (s, y) € t„] x W^. We clearly have that Un-i{s, y) is -measurable and belongs 

to B°°. 

We next consider the difference: 



u{s,y) - Un-i{s,y) = / T{s - r,y - z)a{u{r, z))W{dr,dz) 

^'-ir . (4.4) 

+ I I T{s — r,y — z)b{u{r, z))drdz. 



Jtn-l JR" 

We have the following estimate for the norm in D!"'*^ „ of the above difference: 

Lemma 4.3. For a// m G N and p > 1, there exists a constant C such that, for any s € 
[tn^i,tn], we have: 

\\u{s,y)-un-i{s,y)\t^''' < C - 
Proof. First, observe that Lemma [3^ implies that, for all j = 0, . . . , m: 



sup Et„_, \\D^u{t,z)\\p < C {^{s-tn-i))- , a.s. 

(T,2)G[tn-l,s]xR<* V •^t„_i,sj 

At this point, we can apply Lemma [A. 31 in the following situation: X = u{s,y) — Un-i{s,y), 
^0 = 0, f = b, g = a, V = u, Ii{s,y) = 1, a = in-i> b = s, t = s, j = 0; we point out that 
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here we are doing an abuse of notation, since b in the present section denotes a different type 
of object in comparison with the statement of Lemma I A. 31 (here b is the drift coefficient while 
in the lemma it is simply a time parameter). Hence, one ends up with 

||n(s,y) -'u„_i(s,y)||^-i'^ < C J{s-r)dr^ 

where the constant C depends on {m,p, T) and we recall that J(r) = |^r(r)(^)p/x(dO- 
The latter inequality has been obtained by applying ()3.2p . Therefore, we conclude the proof. 

□ 

Let us start with the decomposition of the term Fn — Fn-i- For this, we set u{X,s,y) := 
Xu{s,y) + (1 — X)un-i{s,y), for A € [0, 1]. Then, by the mean value theorem we have: 



Fn-Fn-i= / T{t - s,x - y)a{un-i{s,y))W{ds,dy) 



+ / T{t - s,x -y)b{un-i{s,y))dyds 

Jt„-i JR'' 

..1 X (4-5) 

T{t- s,x-y){ / a'{u{X,s,y))dX]{u{s,y) -Un-^iis,y))W{ds,dy) 



tn-l 
in 

tn-l 



T{t- s,x-y)[ / 6'(n(A, s, y)) dX ) (n(s, y) - 'u„_i(s, y)) dyds. 



As we will make precise below, the first and second terms on the right-hand side of ()4.5p will be 
called processes of order 1 and 2, respectively, while the third and fourth terms will be called 
residues of order 1 and 2, respectively. 

At this point, we need to introduce some notation, namely we are going to define what we 
understand by processes and residues of order A; € N. The former will be denoted by Jk and 
the latter by Rk- In all the definitions that follow, we assume that s <t. 

For A; = 1, we define: 



Ji{s,t,x)= / T{t - r,x - y)a{un-i{r,y))W{dr,dy). 

J tn-l JR'^ 

For k > 2, the process Jk of order k is defined either of the form: 

s ^ 

Jkis,t,x)= / T{t - r,x - y)a^'^\un-iir,y))T\Jm,{r,r,y)W{dr,dy), 



(4.6) 



where £ < k — 1 and mi + • • • + m£ = A; — 1, or 

Jkis,t,x)= / T{t-r,x-y)b^^\un^i{r,y))Y[jm^ir,r,y)dydr, (4.7) 
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where i < k — 2 and nii + • • • + mi = k — 2; in the case k = 2, the process Jk is defined by 



J2is,t,x)= / T{t-r,x -y)b{un-i{r,y))dydr. 

Jtn-l JK'' 

It is clear that, for any /c € N, the set of processes of order k is finite. Let be an index set 
for the family of processes of order k, so that if a G =24 1 the corresponding process of order k 
indexed by a will be denoted by J^. 

The residues Rk of order k are defined as follows. For k = 1, Ri is defined to be either: 

Ri{s,t,x)= / r{t-r,x-y)( a'{u{X,r,y))dX]{u{r,y) - Un-i{r,y))W{dr,dy) 

Jt„-i JR'' \Jo J 



or 



(4.8) 



Ri{s,t,x)= I T{t-r,x-y)biu{r,y))dydr. (4.9) 

Jtn-l 



For k = 2, set: 



R2{s,t,x)= / T{t-r,x-y)i b' {u{X,r,y)) dX] {u{r,y) - Un-i{r,y))drdy. (4.10) 

Eventually, for any A; > 3, the residue of order k can be one of the following four possibilities. 
First, it can be either: 

Rk{s,t,x) = —1— r f r{t-r,x-y)( f\l-X)'^-'a('^\u{X,r,y))dx) 

[k-iy. Jt„_,jRd \Jo J (4.11) 



X (u(r,y) -Un-i{r,y))''W{dr,dy) 

or 



Rk{s,t,x)= / T{t-r,x-y)a^^\un^i{r,y))T\l^^{r,y)W{dr,dy), (4.12) 



where Imjir,y) is either Rmj{r,r,y) or Jm^{r,r,y), but at least there is one j such that 
Imj{r,y) = h 
can be either: 



Imjir^y) = Rmj{r,r,y). As before £ < k — 1 and mi + ■ ■ ■ + mi = k — 1. Secondly, Rk 



Rk{s, t, x) = r(t - r, X - 2/) V - A)^-26('=-i) (^x(A, r, y)) d>)j 



(4.13) 



X (n(r, y) - ti„_i(r, y))'' ^ dydr 

or 



i?fc(s,i,x)= / / T{t-r,x-y)b^^{un^i{r,y))f\l^^{r,y)dydr, (4.14) 

where £ < k — 2, mi + ■ • • + me = k — 2 and are as in ()4.12p . We will denote here by 
any index set for the residues of order k and R'^ the element corresponding to a G 

Having all these notation in mind, we have the following decomposition for the difference 
Fn — Fn-i- For the statement, let us remind that the drift and diffusion coefficients b and a 
are assumed to be of class 
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Lemma 4.4. For all K (^N we have: 

K 

Fn-Fn-i = Y, Ci{a,k)J^itn,t,x)+ ^ C2 (a, (t„ , t , , (4.15) 

k=l aGsi/^ a^SSii 

where Ci{a^k) and C2{a,K) denote some positive constants. 

Proof. We will use an induction argument with respect to K. First, the case K = 1 follows 
applying the mean value theorem to the function a around the point in the first term on 
the right-hand side of (|4.5p . 

In order to illustrate our argument to tackle the general case, let us prove our statement 
for the case K = 2. In fact, consider the decomposition (|4.5p and observe that it suffices to 
prove that the residue of order 1 (|4.8p can be written as a sum of a process of order 2 and 
residues of order 2. This is achieved by integrating by parts in the integral with respect to dX: 

T{t-s,x-y)[l a\u{X,s,y))dX]{u{s,y) - Un^i{s,y))W{ds,dy) 

tn-l 

tn r 

/ T{t- s,x- y)(j'{un-i{s, y)){u{s, y) - n„_i(s, y))W{ds, dy) 
r / rl 

(2)/ 



+ / / T{t-s,x-y)[ {l-X)a^'^{u{X,s,y))dX\{u{s,y)-Un^i{s,y)YW{ds,dy). 

(4.16) 

Note that, one the one hand, by definition the second term on the right-hand side of ()4.16p is a 
residue of order 2. On the other hand, applying the mean value theorem inside the stochastic 
integral in (|4.4p . one obtains the following decomposition: 

M(s,y) - u„_i(s,?/) = Ji(s,s,y) -h ^ R"{s,s,y). ^^ -^j^^ 

Hence, the first term on the right-hand side of (|4.16p can be written as a sum of a process of 
order 2 and two residues of order 2. This concludes the proof of (|4.15p for K = 2. We point out 
that the same arguments used in the latter final part of the case K = 2 would let us conclude 
that the residue of order 2 (j4.10p can be decomposed as the sum of a process of order 3 and 
residues of order 3. 

Now, we assume that ()4.15p holds for K, and that any residue of order k < K can be 
decomposed as a sum of processes of order k + 1 and residues of order k + 1. Then, we consider 
any residue term of order K in ()4.15p . If this residue of order K is of the form ()4.1ip . by 
integrating by parts one can rewrite it as 



1 



T{t- s,x- y)a^^\u„-i{s, y)){u{s, y) - u„_i(s, y))^W{ds, dy) 



+ j^/(i-s,^-y)(j[ (l-A)V^+i)(tx(A,s,y))dA) 

X (u(s, y) - M„-i(s, y))'^+^W{ds, dy). 

Thus, by ()4.17p . the above expression can be written as a sum of processes of order k + 1 and 
residues of order k + 1. The same computations can be done if the residue of order K is of the 
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form ()4.13p . Eventually, in case the residue is of the form (j4.12p or (j4.14p , we will make use of 
the induction hypothesis. Namely, we will be able to write any Rfn. (K, K, y) therein as a sum 
of processes of order rrij + 1 and residues of order rrij + 1. Therefore, the resulting terms will 
be sums of processes of order k + 1 and residues of order k + 1, which yields that the desired 
decomposition holds for K + 1. This concludes the proof. □ 

— — K 

At this point, we are in position to define the approximation sequence Fn = needed in 
hypothesis (HI). For all K €N, set: 

K 

Fn = An-l{g)^Zn + Fn-l+Y^ ^ Ci{a, k)J^ {tn,t, x) , 

k=i oe^/fc 

where the last term above is the one given in Lemma |4.4[ Note that here we are dealing with 
1-dimensional random vectors, that is random variables. According to the expression of F^ in 
assumption (HI), we have 

K 

Hs,y) ='^{t-s,x-y)a{u„.-i{s,y)) and = ^ ^ Ci(a, A;) J^(t„, t, x), 

k=2 agi4 

and we have taken 7 = ^ • The boundedness of a and the definition of the processes 
guarantee that all the conditions of hypothesis (HI) are satisfied. Let us also remind that 
{Zn, n € N} is an i.i.d. sequence of standard Gaussian random variables which is independent 
of the noise W (see the explanation right after hypothesis (H2d)), and that here we consider 
g{s) = T{t — s), so that 

We next verify condition (H2a). For this, we observe that Lemma 14.41 vields: 
Fn-F„ = -A„_i((7)^Z„+ C2ia,K)R'^{tn,t,x). 

Then, hypothesis (H2a) is a consequence of the following lemma: 

Lemma 4.5. For all m £ N and p > 1, there exists a constant C such that, for any X € N, 
we have: 

||i?i^(t„,t,x)||^;j-i'*" <CA„_i(g)^ a.s. (4.18) 
where Rk denotes any of the four types of residues of order K . 

Proof. We start by assuming that the residue process Rk is of the form ()4.1ip . Then, we will 
apply Lemma Fa. 3 1 in the following setting: = 0, / = 0, g{z) = a^^\Xz + (\ — \)un-i{r,y)), 
Ii{r, y) = u{r, y) — Un-i{r, y) for alH = 1, . . . , io = -f^ and v = u. Note that Lemmas 13.31 and 
14.31 implv that conditions (|A.6p and ()A.7p are satisfied with 7 = i,aj = l,a = Er, respectively 
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(we are again making an abuse of notation since the latter a has nothing to do with the one 
in R'^ above). Altogether, we obtain that: 



\\RKitn,t,x)\t^'''" <C 




< CA„_i(5)'^J\ a.s., 

where we have used the fact that J{r) = f^^ |=^r(r)(^)p/i((i^). 

Let us now assume that the residue Rk is of the form (|4.13p . In this case, we apply 
Lemma 1X3) in the following situation: Xq = 0, g = 0, f{z) = a^^~^\Xz + (1 — A)ti„_i(r, y)), 
Iiir,y) = u{r,y) - M„_i(r,y) for all i = 1, . . . ,io = K - 1, v = u, j = I, Oi = 1, a = K - 1. 
Thus, we end up with: 

K-l 

\\RK{tn,t,x)\t;p'''" <C I f J{t-r)dr] ds 

^t■n. — ^\^t^-, — ^ I 



K-l 



— 1 \ — 1 

<C{tn-tn-l)^n-l{9)^ 

<CA„_i(5()^, a.s. 



In the last inequality we have used the estimate (j3.2p . 

We finally show estimate ()4.18p for residues Rk of the form ()4.12p (the case ()4.14p follows 
along the same lines). This will be deduced thanks to an induction argument on K together 
with an application of Lemma fA.31 We have already checked the validity of (|4.18p for K = 1,2. 
Assume that the residues (j4.12p of order k < K — 1 satisfy the desired estimate. Then, taking 
into account that in ()4.12p there is at least one residue of order k < K hidden in the product 
of Imj{r,y) and Lemma [4. 61 b elow establishes suitable bounds for processes of order K, we can 
conclude by applying again Lemma [A. 31 as follows: Xq = 0, / = 0, g{z) = a^^\z), v = Un-i, 
'y = ^, ai = irij, a = K . Details are left to the reader. □ 

A consequence of the above lemma is that condition (H2a) is satisfied. Moreover, as we 
have already explained at the beginning of the present section, hypothesis (H2b) follows from 
Proposition 14.11 Next we prove (H2c), that is 

Ci<A„_i(5)-W \\his)\\%.ds<C2, (4.19) 

for some positive constants Ci, C2, where we recall that h{s, y) = T{t — s,x — y)a{un-i{s, y)). 
The upper bound is an immediate consequence of [U Theorem 2], since we assume that a 
is bounded. In order to obtain the lower bound in ()4.19p . we can use the same argument 
of [2H Theorem 5.2]. Precisely, let (■i/'e)e be an approximation of the identity and define 
Ge{s,y) := {ipe * h{s)){y), so that Ge{s) € Then, by the non-degeneracy assumption 
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on cr, we can infer that, almost surely: 

/ \\h{s)\\%ds = l\m \\G,{s)f^ds 

Jtn — l Jtn — l 



n-1 
tf-i 



n — 1 



lirn / ds I A{dy) j dz G ^{s , z)G ^{s , z — y) 



>C r I \^T{t-sm\^iJi{di)ds 



'tn-1 

= CA„_i(<7). 

Therefore, we obtain the lower bound in (|4.19|) . which means that (H2c) is satisfied. 
Eventually, we check condition (H2d). To begin with, let us remind that 

K 



k=2 «Gi/fe 

Then, hypothesis (H2d) is proved in the following lemma: 

Lemma 4.6. For all m G N and p > 1, there exists a constant C such that, for any A; € N, 
we have: 

||Jfc(t„,t,2:)||^-i'*" < CA„_i(5)l a.s., 
where denotes an arbitrary process of order k. 

Proof. It is based on several applications of Lemma IA.3[ We shall perform an induction 
argument with respect to k. Assume that Jk is of the form ()4.6p . 

The case k = 1 follows directly from Lemma IA.3I applied to the case Xq = 0, / = 0, 
g{z) = (t{z), V = Un-i and 7 = 0, where the constant on the right-hand side of (|A.6p only 
depends on j,p and T. Similarly, the case k = 2 follows again appealing to Lemma lA.31 with 
-'^o = 0, / = 0, g{z) = cr'(z), v = Un~i, io = I, j = ^, a = 1, and using the result for k = 1. 

Assume that the statement holds for any process of order up to k — 1. Then, if is given 
by (j4.6p . we can apply Lemma [A.3I in the case Xq = 0, / = 0, g{z) = a^^\z), v = Un-i, 
io = -^^ 1 = ^, oii = rui and a = k — 1. Observe that the induction hypothesis guarantees that 
condition (|A.7p is satisfied. Altogether, we obtain: 



\Jk{tn,t, x)\\!i^^p'^" < C 




J{t — s) ds 



<CAn-iig)2, a.s. 



On the other hand, if Jk is of the form ()4.7p . one proceeds similarly as before. Namely, 
we apply Lemma [A. 31 with Xq = 0, g = 0, f{z) = a^^\z), v = n„_i, io = = ^, o.i = f^i, 
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a = k — 2, so that we end up with: 



fc-2 

tn / rs \ 2 



||Jfc(t„,t,x)||^-pi'*" < C / / J{t-r)dr\ ds 

< C {tn - tn-l)An-lig)'^ 

k 

<CA„_i(5()2, a.s., 
where the last inequahty fohows from (|3.2p . This concludes the proof. □ 



With this lemma, we conclude that hypotheses (HI), (H2a)-(H2d) are satisfied and we 
obtain, in view of Definition 12.21 that the random variable F = u{t, x) is uniformly elliptic. 
Therefore, by Theorem 12. 3t we have proved the lower bound in Theorem ll.il 

5 Proof of the upper bound 

This section is devoted to prove the upper bound of Theorem II. li For this, we will follow a 
standard procedure based on the density formula provided by the integration-by-parts formula 
of the Malliavin calculus and the exponential martingale inequality applied to the martingale 
part of our random variable u{t, x) for (t, x) G (0, T] x M'^ (this method has been used for 
instance in |18l Proposition 2.1.3] and [3 [5]). We remind that we are assuming that the 
coefficients b and a belong to ^{,°°(M) and the spectral measure /i satisfies 



Moreover, we have that: 



u{t,x) = Fo+ [ I T{t- s,x-y)cj{u{s,y))W{ds,dy) 

Jo JR<* 



+ 



/ / ^{t- s,x-y)b{u{s,y))dyds, a.s., 
Jo Jr'' 



(5.1) 



where Fq = (r(i) * uo){x). 

To begin with, we consider the continuous one-parameter martingale {Ma,c^a, < a < t} 
defined by 

Ma= [ [ T{t-s,x-y)a{u{s,y))W{ds,dy), 



Jo JW^ 

where the filtration {=^a, < a < i} is the one generated by W. Notice that Mq = and one 
has that 

{M)t = \\T{t - ■,x - *)a{u{-,*))\\.M- 

Since a is bounded, we clearly get that {M)t < C2^{t), a.s. for some positive constant C2 (see 
for instance [H Theorem 2]). 

On the other hand, since the drift b is also assumed to be bounded and r(s) defines a 
probability density, we can directly estimate the drift term in (|5.ip as follows: 







T{t - s,x- y)b{u{s, y)) dyds 



< C3 T, a.s. (5.2) 
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We next consider the expression for the density of a non-degenerate random variable that 
follows from the integration-by-parts formula of the Mallavin calculus. Precisely, we apply [181 
Proposition 2.1.1] so that we end up with the following expression for the density pt^x of x): 



Pt,x{y) = E '^{u{t,x)>y}^{Du{t,x) \\Du{t,x)\\J^ 



y e 



where 5 denotes the divergence operator or Skorohod integral, that is the adjoint of the Malli- 
avin derivative operator (see [181 Ch. 1]). Taking into account that the Skorohod integral 
above has mean zero, one can also check that: 

Pt,x{y) = ^{u{t,x)<y}S{Du{t,x)\\Du{t,x)\\^^) , yeR. 

Then, owing to (jS.ip . Proposition 2.1.2] and the estimate (|5.2p . we can infer that: 

Pt,xiy) < Ca,/3,gP{|Mt| > \y - Fo\ - CgT}!/" 

X (E[\\Du{t,x)\\^] + \\D^u{t,x)\\^^(^^.j^^^2-^\\\\Du{t^ ^ ^ 

where a,l3,q are any positive real numbers satisfying ^ + ;g + | = 1- Thus, we proceed to 
bound all the terms on the right-hand side of ()5.3p . 

First, by the exponential martingale inequality (see for instance [181 Sec. A2]) and the fact 
that (M)t < C2$(t), we obtain: 

P{lAf,| >\y- Fol - C2T} < 2exp(^-fc^^^^) . (5.4) 

Secondly, we observe that the following estimate is satisfied: for all p > 0, there exists a 
constant C, depending also on T, such that 

E{\\Du{t,x)\\J^)<Cm-''- (5.5) 

Remark 5.1. Indeed, this bound could be considered as a kind of particular case of ()4.ip . 
Nevertheless, though its proof is very similar to that of Proposition 14. II and we will omit it (an 
even more similar proof is given in [21' Theorem 6.2]), it is important to make the following 
observation. Namely, the fact that in ()5.5p . compared to Proposition l4.1l we are not considering 
a time interval of the form [t„_i,f„] but directly [0,t], makes it possible to obtain the desired 
estimate ()5.5p without using condition (|3.2p . Therefore, owing to Remark 14.21 one could verify 
that the density upper bound we are proving in the present section turns out to be valid for 
the stochastic wave equation in space dimension d € {1, 2, 3}, with the only differences that Fq 
and the term C3T should be replaced by the corresponding contribution of the initial conditions 
and csT^, respectively. Note that the power in that case is due to the total mass of the 
associated fundamental solution considered as a measure in (see [H Example 6]). Finally, 
we point out that, for spatial dimension one, density upper bounds for a reduced stochastic 
wave equation have been obtained in [6]. 

Going back to our estimate ()5.5p . let in particular p = \ and p = (3 m. such a way that, 
respectively: 

E[||Z)7x(t,x)||,^^] <C<i>(t)-i/2 and \\\\Du{t,x)\\J^\\L,^^^ < C ■ (5-6) 
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Eventually, Lemma 13.31 implies that 

||Z)\(t,x)||^,(^.^»2) < C^{t) < (5.7) 

where the latter constant C depends on T. Hence, plugging estimates (|5.4p - (j5.7p into ()5.3p we 
end up with: 

.*(*)-'^exp(-i-||^), 

where the constants Cj do not depend on {t,x). This concludes the proof of Theorem ll.il 

A Appendix 

The first part of this section is devoted to recall the construction of the Hilbert-space- valued 
stochastic and pathwise integrals used throughout the paper, as well as establish the corre- 
sponding conditional L^-bounds for them. In the second part, we state and proof a technical 
result that has been very useful in the proofs of Section HI 

A.l Hilbert-space- valued stochastic and pathwise integrals 

Let us briefly explain the construction of stochastic and pathwise integrals in a Hilbert-space- 
valued setting, the former of course being with respect to W. This is an important point in 
order to consider the linear stochastic equation satisfied by the iterated Malliavin derivative 
of the solution of many SPDEs (for a more detailed exposition, see \26\ Section 2] and [28]). 

More precisely, let £/ he a separable Hilbert space and {K{t,x), {t,x) € [0, T] x W^} an 
=«2/-valued predictable process satisfying the following condition: 

sup ii;(||i^(t,x)||^)<+oo, (A.l) 
where p > 2. We aim to define the ^/-valued stochastic integral 

G-Wt= f [ G{s,y)W{ds,dy), t e [0,T], 

Jo JR-i 

for integrands of the form G = T(s,dy)K(s,y), where we assume here that T is as described 
in Remark 13.41 In particular, T satisfies condition ()1.4p . that is: 

r [ \,^Tit){0\'^fi{dOdt < +00. 

Jo JlRd 

Note that these assumptions imply that G is a well-defined element of L'^{n x [0, T]; (g) i^). 
Recall that we denote by {^t, t > 0} the (completed) filtration generated by W. Then, 
the stochastic integral of G with respect to W can be defined componentwise, as follows: let 
{ej, j S N} be a complete orthonormal basis of ^ and set G^ := V{s,dy)K^ {s,y), where 
K^{s,y) := {K{s,y),e,)^, j G N. We define 

G-Wt:=Y,G^ -Wt, 

j&i 
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where - Wt = Jq J^d r(s, y)K^{s, y)W{ds, dy) is a well-defined real- valued stochastic integral 
(see [261 Remark 1]). By (|A.ip . one proves that the above series is convergent in and 
the limit does not depend on the orthonormal basis. Moreover, {G ■ Wt,.^t, t € [0, T]} is a 
continuous square-integrable martingale such that 

E{\\G - WtWI,) = E{\\G\\%.^^,^). 

We also have the following estimate for the pih. moment of G ■ Wt (see |26i Theorem 1]): for 
aU t € [0,r], 

E{\\G-Wtr^)<Cpm'^'^ I sn^E{\\K{s,x)rjJ[s)ds, (A.2) 
where we remind that 

m= f I \^T{sm\'' iiidOds and J{s)=f |^r(.)(e)| ^(^0- 

Next, we consider a conditional version of (|A.2p : 

Lemma A.l. For all p >2 and < a < b < T , we have: 

l-b 

E[\\G-Wi,-G-Wa\U^a]<Cp{^b)-^a))^-^ / snv E[\\K{s,x)fJ^a]J{s)ds, a.s. 

The proof of this result is essentially the same as its non-conditioned counterpart (|A.2p . 
except of the use of a conditional Burkholder-Davis-Gundy type inequality for Hilbert-space- 
valued martingales. 

Let us now recall how we define the Hilbert-space- valued pathwise integrals involved in the 
stochastic equations satisfied by the Malliavin derivative of the solution. Namely, as before, 
we consider a Hilbert space i/, a complete orthonormal system {e^, j G N}, and an .^/-valued 
stochastic process {Y{t,x), {t,x) G [0,T] x R'^} such that, for p>2, 

sup Ei\\Y{t,x)r^)<+^. (A.3) 

Then, we define the following pathwise integral, with values in L'^{Q] jz/): 

A:= I I Y{s,y)T{s,dy)ds ■.= y^( [ [ {Y{s,y),ej)^T{s,dy)ds] ej, t G [0,T], 
Jo Jr-i \Jo jr<* / 

where F is again as general as described in Remark 13.41 Moreover, a direct consequence of the 
considerations in [251 p. 24] is that: 

E{\\A\C)<( tns,^'')dsY [\npE{\\Y{s,z)rjr{s,R'')ds. (A.4) 

\J0 J Jo zGRrf 

In the paper, we need the following straightforward conditional version of the above estimate 
Lemma A.2. Let p >2. Then, for any a-field 'S , we have: 

E[\WW^]< ( fris,R'')dsY f snpE[\\Yis,z)rj';^]Tis,R^)ds, a.s. 
\Jo J Jo zeM'* 
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A. 2 An auxiliary result 

Let a < b <t and x € M.'^, and consider the following random variable: 

X = Xo{t,a,b)+ / T{t-s,x-y)fiv{s,y))T[l,{s,y)dyds 



r-b r- «0 

+ / / T{t - s,x - y)g{v{s,y))\\h{s,y)W{ds,dy), 



(A.5) 



where Xo(t,a, 5) is a ^a-measurable random variable, /, (7 : M — )• M are smooth functions and 
V and li are certain smooth stochastic processes (in the Malliavin sense). As usual, F denotes 
the fundamental solution of the stochastic heat equation in M*^. In the next result, we provide 
an estimate for the pth moment of the iterated Malliavin derivative of order m for the random 
variable X. 

Lemma A. 3. Assume that f,g & '^^°°(M) and that v and are smooth stochastic 

processes on [a,b] for which the random variable (jA.Sp is well-defined. Let m € N and p > 2. 
Suppose that, for all j = 0, ...,m and s E [o-jt], there exist constants c{j,p,a,s) > 1 and 
C{m,p) such that the former increases in p and 

Ea\\D^v{s,y)f^^^^ <c{j,p,a,s) and (A.6) 
EjDUiis,yW^y<Cirn,p)(rJit-r)dry \ (A.7) 



for all i = l,...,io and y € M"', and some 7 > 0, ai, . . . ,aig > 0, where we recall that 
J{t-r) = J^,\^r{t-rm\^t,{d^). 

Then X G J}"^'°° and the following is satisfied: if we set q := ai + • • • + ai^, there exist p' 
and p* such that 



1 1 
PI „ \ p 



EaWD^^Xf^^^^^j < ^EjD^Xo{t,a,b)\,^^ 

/b / /-s \ 07 

[c*{m-l,p',a,s)'^' + c{m,p*,a,s)] i J{t-r)dr\ ds 

+ 02(^1^ J{t-r) dry [c*{m-l,p\a,sr + c{m,p\a,s)] 



(A.8) 



a'yp 



J{t - r) drj J{t - s) ds 

where Ci,C2 are some positive constants possibly depending on m and p, such that if f = 
then Ci = 0, and if g = then C2 = 0. We also use the notation c*{m — l,p',a,s) := 
maxo<j<m-i c(i)P') s) and we set c{—l,p,a, s) =0. 
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In the case where c{j,p' , a,s) < 1 for all j = 1, . . . , m, estimate (jA.SP is replaced by 



<{Ea\\D"^Xo{t,a,b)r^, 



+ Ci [c* {m — l,p' ,a, s) + c{m,p* ,a, s)] ( / J{t — r)dr 



a7 



ds 



+ C2 



J{t — r) dr 



2 p 



[c*(m — a, s) + c{m,p* , a, s)]' 



J{t -r)dr] J{t - s) ds 



Proof. The assumptions on the functions /, g and the processes v, li clearly yield that X € 
pm,oo_ fjgj^gg^ proof will be devoted to establish estimate (jA.Sp . The method is similar to 
that of Lemma 14 in |10j . and thus we will only focus on the parts of the proof which really 
exhibit a different methodology. 
To begin with, we have that 



Ea\\D'^X\f, 



< [Ea\\D"'Xo{t,a, 



+ / T{t-s,x-y)[Ea\\D'^f{v{s,x-y))\\h{s,x- 

T{t- s,x- y)g{v{s, y)) Ii{s, y)W{ds, dy) 



dyds 



(A.9) 



+ Ea 



i=l 



The second term on the right-hand side of the above inequality may be bounded in the same 
manner as for the corresponding term in the proof of [10| Lemma 14] (see (13) therein). 
Indeed, one makes use of the Leibniz rule for the iterated Malliavin derivative and applies a 
generalization of Holder's inequality, which altogether yields to 



/ / V(t-s,x-y)\Ea\\B^f{v{s,x-y))T\h{s,x 



a.b 



dyds 



(A.IO) 



<C{m,p) / [c* {m — l,p ,a, s)^ + c{m,p* ,a, s)) i / J{t — r)dr\ ds. 



a-y 



Here, we have p' := max^p^, where this maximum is taken over a finite set whose cardinal 
depends on m, and pk > denote a certain set of real numbers needed for the application 
of the above-mentioned Holder's inequality (see ^Oi p. 455] for details). On the other hand, 
p* := p in the case where /j is constant, for all i; otherwise p* := p' . Furthermore, in the case 
where c{j,p' , a, s) < 1 for all j = 1, . . . , m, estimate (jA.lOp turns out to be 

/b / /-s \ "7 

[c*{m — l,p ,a,s) + c{m,p* ,a,s)) i / J{t — r)dr\ ds. 
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In order to deal with the last term on the right-hand side of ()A.9p . the computations slightly 
differ from those used in Lemma 14], since we are considering the more general setting 
determined by the Hilbert space J^. We will use the following notation: s := (si, . . . , Sm) G 
[a, b]'", and := (si, . . . , Sj+i, . . . , Sm)- Then, for instance, for any smooth random 
variable Y, D^Y denotes the ^'^'"-valued random variable defined by (L>™y)(s,*). 

Using this notation, we can first infer that, for all s £ [a, b]"^: 



T{t - s,x- y)g{v{s, y)) JJ Ii{s, y)W{ds, dy) 



j=0 



+ 



f [ Tit-s,x-y)D^(givis,y))f\li{s,y))w{ds,dy), 



where these equalities are understood as random variables with values in J^'^'^, and we recall 
that -k denotes the J^- variable. Thus, the last term in (|A.9|1 may be bounded by Ai +A2, with 



Al = Ea 



{a,by' 



^0 



b 



\ i=0 

r{t -s,x- y)Df (^g{v{s, y)) f[ I,{s, y) j W{ds, dy) 



aVsiV---VSm 





V 
2 




ds 










2 


W{ds,dy) 


ds 









The term Ai can be treated using similar arguments as those of |27^ Lemma 2]. Indeed, if we 
set 



Zj{s,y) :-- 



D'^^Mg{v{s,,y))\{h{s,,y) 



i=0 



_^®(m-l) 



{a,by 



(which is a real- valued random variable), then the fact that F is a smooth function and Cauchy- 
Schwarz and Holder's inequalities yield, up to some positive constant, 

A^ < Ea / dsj / A{dy) / dz T{t — Sj,x — z)T{t — Sj,x — z + y) 
Zj{s, z)Zj{s,z - y) ds{j) 

' f-b \ "I i-b / /• \ 2 

/ J{t-r)dr) V / sup Ea / \Zj{s, y)\^ dsiJ) J(t - sj) dsj. 

Ja J J a yeRd yj(a,6)™-l / 

(A.ll) 

The §th moment appearing in the latter term can be estimated using the same arguments we 
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have commented above to obtain (jA.lOp . Namely 



Ea 



\Zj{s,y)\^ ds{j) 



= eA D"'-^[g{v{s,,y))\\h{s,,y) 
,'ia,br-' J \ \ f=l . 

< C{m-l,p)c*{m-l,p',a,Sj)^"'-^'^P ( [ ' J{t-r)dr 
Plugging this bound in (|A.lip we eventually end up with 

< mC{m — l,p) J{t — r) dr 

X / c*(m- l,p',a,s)(™-i)f / J{t-r)dr] J{t-s)ds 



5(m-l) 



a, 6 
«7P 



£-1 
2 ^ 



(A.12) 



In order to bound the term l^t us apply Lemma lA.ll in the particular case where 



■^T and 



G = D™h7(^;(s,y))n^i(s,2/) 



Hence, we have that, up to some positive constant. 



J{t — r) dr 



sup Ea 



D^{g{v{s,y))\{h{s,y) 



J{t — s) ds. 



At this point, one applies the same method that we have used to obtain estimate (jA.lOp (see 
also the last part in the proof of |10^ Lemma 14]), so that we can infer that 



b \ f-l rb 



J{t - r) dr 



[c*(m — l,p', a, s)™ + c{m,p* ,a, s)]' 



J{t - s) ds. 



We conclude the proof by putting together estimates (|A.10p , (|A.12p and (|A.13|1 . 



(A.13) 



□ 



Remark A. 4. As for Lemma 13.31 the above lemma still remains valid for a slightly more 
general situation. Namely, in the case where T satisfies the assumptions specified in Remark 
3.41 such as for the stochastic wave equation in space dimension d G {1, 2, 3}. In such a general 
setting, the proof of Lemma I A. 3 1 becomes even more technical and tedious since one needs to 
smooth r by means of an approximation of the identity. For the sake of conciseness, we have 
decided to focus the proof on our stochastic heat equation. 
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